Accurate retrieval of structural information from laser-induced photoelectron and 
high-harmonic spectra by few-cycle laser pulses 
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By analyzing "exact" theoretical results from solving the time-dependent Schrodinger equation of 
atoms in few-cycle laser pulses, we established the general conclusion that differential elastic scat- 
tering and photo-recombination cross sections of the target ion with free electrons can be extracted 
accurately from laser-generated high-energy electron momentum spectra and high-order harmonic 
spectra, respectively. Since both electron scattering and photoionization (the inverse of photo- 
recombination) are the conventional means for interrogating the structure of atoms and molecules, 
this result shows that existing few-cycle infrared lasers can be implemented for ultrafast imaging of 
transient molecules with temporal resolution of a few femtoseconds. 
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Electron diffraction and X-ray diffraction are the con- 
ventional methods for imaging molecules to achieve spa- 
tial resolution of better than sub- Angstroms, but they are 
incapable of achieving temporal resolutions of femto- to 
tens of femtoseconds, in order to follow chemical and bio- 
logical transformations. To image such transient events, 
large facilities like ultrafast electron diffraction method 
[H and X-ray free-electron lasers (XFELs) are being de- 
veloped. Instead of pursuing these expensive technolo- 
gies, here we provide the needed quantitative analysis to 
show that existing few-cycle infrared lasers can be imple- 
mented for ultrafast imaging of transient molecules. 

When an atom is exposed to an infrared laser, the atom 
is first tunnel ionized with the release of an electron. This 
electron is placed in the oscillating electric field of the 
laser and may be driven back to revisit its parent ion. 
This re-encounter incurs, in the second step, various elas- 
tic and inelastic electron-ion collision phenomena where 
the structural information of the target is embedded. The 
possibility of using such laser-induced returning electrons 
for self-imaging molecules has been discussed frequently 
in the past. Theoretical studies of laser-induced elec- 
tron momentum images of simple molecules do show in- 
terference maxima and minima typical of diffraction im- 
ages, but they are observed only for large internuclear 
distances [2, 3, 4, 5, ii]. In addition, for quantitative anal- 
ysis, the role of laser fields on these diffraction images still 
has to be understood . More recently, it was reported 
that the outermost molecular orbital of N2 molecules can 
be extracted from the high-harmonic generation (HHG) 
spectra using the tomographic procedure This inter- 
esting result has generated a lot of excitement, but the 
reported results are obtained based on a number of as- 
sumptions [^, 0, [lo| . To make dynamic chemical imaging 
with infrared lasers as a practical tool, general theoretical 
considerations, especially the validity of the extraction 
procedure, should be examined carefully. 



In this Letter, we show that elastic scattering cross 
sections of the target ion by free electrons can be ac- 
curately extracted from laser-induced photoelectron mo- 
mentum spectra. We also show that accurate photo- 
recombination cross sections of the target ion can be 
extracted from the HHG spectra. Our conclusions 
are based on "exact" theoretical results by solving the 
time-dependent Schrodinger equations (TDSE) of sev- 
eral atoms in intense laser fields. While our conclusions 
are derived from atomic targets, the same conclusions 
should hold for molecules too (where "exact" calcula- 
tions are not possible). For molecules, the results have 
important implications. Both elastic scattering and pho- 
toionization are the standard methods for studying the 
structure of atoms and molecules in conventional energy 
domain measurements, thus high-energy photoelectrons 
and high harmonics generated by infrared lasers offer the 
promise for revealing the structure of the target, with the 
temporal resolution offered by the ultrashort laser pulses. 

Consider a typical few-cycle laser pulse, with mean 
wavelength of 800 nm and peak intensity of 10^"* W/cm^. 
The electric field F{t) — —dA{t)/dt and the vector po- 
tential A{t) of such a laser pulse are depicted in Fig. 
1(a). By placing a hydrogen atom in such a laser pulse, 
we solved the TDSE to obtain the photoelectron energy 
and momentum distributions, shown in Figs. 1(b) and 
1(c), respectively. Fig. 1(d) shows the electron momen- 
tum image of Ar in the same laser pulse. The theoretical 
method for solving the TDSE has been described previ- 
ously [11 1 . 



In Fig. 1(b), two particular energies, 2Up and IQUp, are 
marked, where [/p = A§/4 is the ponderomotive energy, 
with Aq being the peak value of the vector potential of 
the laser pulse. (We use atomic units in this paper.) 
These two energy values are important, the former is the 
maximum energy the electron can reach if it is released 
by the laser field alone, while the latter is the maximum 
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FIG. 1: (Color online) (a) Schematic of the electric field (E) 
and the vector potential (A) of a typical few-cycle pulse, (b) 
Energy spectra of atomic hydrogen ionized by a 5 fs (FWHM) 
laser pulse, with mean wavelength of 800 nm and peak inten- 
sity at 10^'' W/cm'^. (c) Normalized 2D photoelectron mo- 
mentum spectra of atomic hydrogen. The images are renor- 
malized for each photoelectron energy to reveal the global 
angular distributions, (d). Momentum images of Ar in the 
same pulse. 



energy the electron will have if it is back scattered by the 
parent ion [l^]. To display the full electron momentum 
image surface in a single plot, in Figs. 1(c) and 1(d) we 
normalized electron momentum distributions such that 
the total ionization yield at each electron energy is the 
same. We have chosen the horizontal axis to be along the 
direction of the laser's polarization and the vertical axis 
along any direction perpendicular to it (due to cylindrical 
symmetry of the linearly polarized hght). 

In Figs. 1(c) and 1(d), we note that in the large mo- 
mentum region, each image exhibits two half circular 
rings, one on the "left" and one on the "right" , with the 
center of each circle shifted from the origin. The rings are 
very similar for the two targets. We will call these cir- 
cular rings back rescattered ridges (BRR), representing 
electrons that have been rescattered into the backward 
directions by the target ion. The BRR on the "right" is 
from electrons born at time near 'a' (Fig. 1(a)), travel- 
ing to the "right" and then returning to the target ion 
at time near 'b', where they are rescattered back to the 
right. Each momentum half circle is represented approxi- 
mately by ArPz +PoPr, where the second term is the mo- 
mentum of the backscattered electron and the first term 
is the momentum added to the electron as it propagates 
from 'b' to the end of the laser pulse. The magnitude 
of the momentum pq is related to the ponderomotive en- 
ergy by 3.17C7p = pI/2 (where Up = Aj/A, and Ar is the 
vector potential at 'b' ), which is the maximum energy of 
electrons that return to revisit the parent ion. For back 
scattered electrons, the two momentum terms add to give 
high-energy photoelectrons, reaching a maximum of lOC/p 
for electrons that have been scattered by 180° 12]. If the 
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FIG. 2: (Color online) Angular distributions of photoelectrons 
along the BRR compared to the differential elastic scattering 
cross sections of the target ion. Each BRR is taken to be 
the outer half circle on the right side of Fig. 1(c) and 1(d), 
respectively, (a) for H target, (b) for Ar, (c) for Ne, (d) for 
Xe. 

electrons are scattered into the forward direction, the two 
momentum terms subtract from each other, resulting in 
lower energy electrons. Similar back scattered electrons 
are found on the "left". These are from electrons that 
were born near 'a" and rescattered back to the left near 
'b" [Fig. 1(a)]. Both the shift of the center and the radius 
are smaller due to rescattering occurring in the smaller 
vector potential near 'b". 

In Figs. 1(c) and 1(d), we note that the yields on the 
BRR for H and Ar are quite different, where the former 
is monotonic and the latter has a clear minimum. Taking 
the actually calculated photoelectron yields (without the 
normalization as shown in the figure) we compare the an- 
gular dependence of the intensities along BRR with the 
elastic differential cross sections of the target ion by free 
electrons at energy E = Pq/2. The results are shown 
in Fig. 2(a) for H target and 2(b) for Ar, where the 
scattering angles are measured from the direction of the 
"incident" electron beam. Good agreement between the 
two results for each target atom can be seen. Such good 
agreement has been duplicated at different laser intensi- 
ties and other atomic systems [see Ne and Xe in Figs. 2(c) 
and 2(d), respectively.]. These results show that laser- 
induced momentum images on the BRR can be used to 
obtain elastic scattering cross sections of free electrons by 
the target ion. For atomic hydrogen, the elastic scatter- 
ing cross section is given by the Rutherford formula. For 
other atomic ions, minima in the differential cross sec- 
tion do occur and they are due to interference between 
electrons scattered by the short-range potential and by 
the Coulomb potential [13]. We comment that momen- 
tum spectra due to back scattered electrons have been 
observed experimentally earlier (14| but there has been 
no quantitative analysis. 

According to the intuitive rescattering model, the pho- 
toelectron yield along the BRR may be interpreted as due 
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FIG. 3: (Color online) Electron wave packets extracted from 
photoelectron momentum images along the BRR are shown to 
be identical for scattering angles of 160° and 180°, extracted 
from (a) H target, (b) Ar target, using the same ionizing laser, 
(c) The extracted electron wave packet from H target is the 
same over the angular range of 130° to 180°. 



to back scattering of the returning electron wave packet. 
To test this idea, we write the photoelectron momentum 
yields I{p) along the BRR by I{p) = a{po,0)F{pa,e), 
with p — ArPz +poPr, where a{po,9) is the elastic dif- 
ferential cross section for each ion by a free electron with 
energy E — Pq/2, and 9 is the scattering angle of the free 
electron. 

In Fig. 3(a), we show that the extracted F{po, 9) from 
H at two different angles. They are essentially identical 
such that we may identify F{pq) = F{po, 9) as the wave 
packet of the returning electrons with momentum near 
Pq. Note that the wave packet extracted from Ar target, 
as shown in Fig. 3(b), is essentially identical to Fig. 3(a) 
except for a small shift of the center from po = 1.22 
to 1.25. The width of the wave packet is found to be 
independent of the target. In Fig. 3(c) we show that 
P{Po) — F{po,9) indeed holds well for electrons that 
have been back scattered for angles larger than about 
130°. Note that a separate electron wave packet can 
be retrieved from the photoelectron momentum spectra 
measured on the "left" . The wave packet analyzed above 
is for a 5 fs (FWHM) pulse with the carrier-envelope 
phase (CEP) (j) = 0. The returning electron wave packet 
depends on the CEP and can be used to measure the CEP 
of a few-cycle pulse. For long pulses, electrons along each 
side of the BRR will exhibit oscillations characteristic of 
ATI peaks, and elastic scattering cross sections by free 
electrons can be extracted from the envelope of the mo- 
mentum images on the BRR [15]. 

The above results, based on the exact solution of 
TDSE, clearly established that the electron yields on the 
BRR, can be viewed as the backscattering of the return- 
ing electron wave packet. According to the "simple man" 
model, the HHG is due to the emission of photons follow- 
ing the photo-recombination of the same returning elec- 



trons. Since the returning electron wave packets have 
been shown largely independent of target atoms for the 
same laser pulse, the difference in the HHG spectra would 
then be attributed to the recombination cross sections. 
To check this idea, we choose a companion atom to have 
the same ionization potential Ip such that the cutoff ener- 
gies will be identical. In the calculation, we compare the 
HHG spectra for Ne and scaled model hydrogen atom, 
generated by the same laser field, by solving the TDSE. 
For Ne atom, it is ionized from the 2po state, and for H, 
we chose an effective charge such that its ground Is state 
has the same ionization potential as Ne. For the laser, 
we chose a 5.2 fs pulse, with mean wavelength of 1064 
nm and intensity of 2 x 10^^ W/cm^. The longer wave- 
length was used such that harmonics would extend over 
a larger photon energy range. From the calculated HHG 
yield, we divide each by its "exact" photo-recombination 
cross section. We compare in Fig. 4(a) the resulting 
"electron wave packets" (normalized) vs the HHG or- 
der. Clearly the two electron wave packets are very close 
to each other, showing that the differences in the HHG 
spectra between the two targets are coming from the dif- 
ferent photo-recombination cross sections. Note that the 
kinetic energy of the photoelectron is related to the pho- 
ton energy by hcu = + Ip. Also note that the de- 
duced "electron wave packet" shows many oscillations. 
To study the structural information of a target, it is much 
easier to factor out these oscillatory features by using a 
known companion atomic target. 

The good agreement in the deduced "electron wave 
packet" in Fig. 4(a) prompts us to ask whether one can 
obtain accurate photo-recombination cross sections from 
the HHG spectra of an unknown target by comparing it 
with the HHG spectra from a known atom with identi- 
cal or nearly identical ionization potential. For our pur- 
pose we generated the HHG spectra for Ar and scaled 
H atoms, with the effective charge of H chosen to have 
identical ionization potential as Ar. Different laser inten- 
sities and wavelengths were used to generate HHG from 
which Ar recombination cross sections are derived. If 
the procedure is valid, the extracted recombination cross 
sections should be independent of the lasers used, except 
for the range of photon energies covered. In Fig. 4(b) 
we compare the extracted cross sections from the HHG 
generated by various laser pulses with the "exact" photo- 
recombination cross section calculated for Ar. One can 
see that the deduced values scatter nicely around the "ex- 
act" one. The fluctuation of the extracted cross sections 
can be reduced if the HHG intensity is taken from the 
averaged HHG amplitudes with nearby laser intensities. 
Thus the smooth black line is obtained from averaging 
over eleven intensities within ±5% of the mean intensity 
of 2 X 10^'' W/cm^. This coherent averaging sharpens the 
odd harmonics and reduces the harmonic yields in be- 
tween, similar to the effect of propagation of HHG in the 
medium. Note that in the "exact" calculation we include 
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FIG. 4; (Color online) (a) Comparison of the "electron 
wavepackets" extracted from the HHG spectra of Ne and 
scaled H generated by a 5.2 fs laser pulse with peak inten- 
sity of 2 X 10^** W /err? and mean wavelength of 1064 nm. 
Extracted photo-recombination cross sections of Ar (b), Ne 
(c) and Xe (d), using different laser pulses, plotted vs the 
photon energy. For each panel, the dashed red line is the 
"exact" theoretical cross section, calculated using the same 
model potential, as used in the TDSE. (Jo = 10^"* W/cm^.) 

dipole transitions from 3p to both s and d continuum 
states. The minimum in the cross section occurs near 
the Cooper minimum where the d-wave electric dipole 
moment changes sign. Comparing the two curves we can 
say that accurate photo-recombination cross sections in- 
deed can be extracted from the HHG yields. Similar test 
showing good agreement has been made also on Ne and 
Xe atoms (Figs. 4(c) and 4(d)), again using scaled hy- 
drogen as the companion atoms. 

Following last paragraph, we now comment the as- 
sumptions made in Itatani et al Q- Their assumption 
that the returning electron wave packets for Ar and N2 
(they have nearly identical ionization potentials) are the 
same under the same laser pulse is confirmed by our cal- 
culations which are based on exact solution of the TDSE. 
However, it was also assumed 0, that the dipole ma- 
trix elements can be calculated using plane waves instead 
of "exact" scattering waves. This leads to results that are 
quite different from ours. 

In this Letter we have identified the spectral region 
where the nonlinear laser-atom interaction can be simpli- 
fied to extract the linear interaction between a returning 
electron wave packet with the atomic ion. Even though 
the results were presented for atomic targets, the same 



simplifications should hold for molecular targets as well. 
For molecules, this opens up the exciting opportunity of 
using infrared lasers for ultrafast imaging of molecules 
undergoing structural transformation. Both elastic elec- 
tron scattering and photoionization (the inverse of photo- 
recombination) processes are the well-tested means for 
probing the structure of molecules. They depend lin- 
early on the electron current and photon intensity, re- 
spectively. Theoretically, these cross sections can be cal- 
culated very precisely. A computer package for obtaining 
elastic scattering cross section and photoionization cross 
section within such a model has been published recently 
by Tonzani By extracting elastic and/or photo- 

recombination cross sections of molecules using few-cycle 
infrared lasers, structural changes of the molecules can 
be determined with temporal resolution of a few fem- 
toseconds. In conclusion, we have established the theo- 
retical foundation for carrying out structural analysis of 
molecules with infrared lasers. If this roadmap is imple- 
mented experimentally, table-top infrared lasers would 
offer a very competitive new technology for ultrafast 
time-resolved chemical imaging, with temporal resolution 
down to a few femtoseconds. 
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